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SUMMARY 

The small -perturbation equations of motion of a flexible aircraft with an active control 
technology (ACT) system were developed to evaluate the stability and performance of the 
controlled aircraft. The total aircraft system was formulated in state vector format and the 
system of equations was completed with fully unsteady and low-frequency aerodynamics for 
arbitrary, complex configurations based on a potential aerodynamic method. The ACT sys- 
tem equations have been incorporated in the digital computer program FCAP (Flight Control 
Analysis Program) which can be used for the analysis of complete aircraft configurations, 
including control system, with either low-frequency or fully unsteady aerodynamics. The 
application of classical performance analyses including frequency response, poles and zeros, 
mean-square response, and time response in FCAP in state vector format was discussed. 

INTRODUCTION 

The integrated study of the interactive effects of the flight control system in the active 
control of flexible aircraft has received considerable attention in recent years. In particular. 
Active Control Technology (ACT) is being investigated for improving ride quality, decreas- 
ing structural deformation, extending the fatigue life of the aircraft, relaxing static stability 
requirements, suppressing flutter, and reducing structural loads. 

A new computer program. Flight Control Analysis Program (FCAP) has been developed 
for NASA to analyze ACT systems (refs. 1 and 2). The program was designed in a modular 
fashion to incorporate aircraft dynamics, aerodynamics for complex configurations, and 
sensor, actuator, and control logic dynamics, os well as analysis methods for determining 
stability and performance of the ACT system. The formulation of the total aircraft dynamic 
system for FCAP was unified by casting all the equations in state space format. This paper 
presents the state-vector formulation of the ACT system, and discusses its application in 
FCAP for the performance analysis of ACT systems. 


*This paper was derived from work conducted under NASA Contract NAS 1-13371 . 

**Consultant to Aerospace Systems, Inc. 
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Subscripts: 
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state-space matrices of ACT system and subsystems 

matrix defined by Equation (17) 

matrix defined by Equation (16) 

matrix of transfer functions 

imaginary part of complex frequency, s 

generalized mass/inertia matrix 

dynamic pressure 

Lagrangian generalized coordinates , 
covariance matrix of ACT system outputs 
output vector for ACT system and subsystems 
complex frequency 

covariance matrix of ACT system inputs 
generalized aerodynamic force coefficients 

coefficient matrices of aircraft dynamics equations at sensor locations 

input vector for ACT system and subsystems 

matrix of pilot and guidance system commands 

generalized aerodynamic forces in uniform flow 

matrix of aerodynamic forces due to turbulence 

covariance matrix of ACT system state variables 

state vector 

actuator 

aircraft displacement 
control logic 
aircraft rate 
sensor 
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Superscripts; 

T transpose 

(0), (1) coefficients of power series expansion 

A tilda (~) over a symbol indicates that it is designated in the Laplace domain. A dot 
over a variable indicates time differentiation. 

ACT ANALYSIS 

The ACT system formulated in FCAP is shown in Figure 1 . External disturbances to the 
ACT system are seen to be atmospheric turbulence and gusts contributing to the aerodynamic 
forces and moments, and pilot or guidance system commands introduced through the control 
logic. The aircraft dynamic system includes both rigid-body and flexible-body dynamics. 

The analysis of ACT systems is unified by casting all system equations in either the time 
domain or the frequency domain, and in similar format. In state space methods, the motion 
of a given dynamic system is described by the following pair of matrix equations (ref. 3): 

X = A X + B u 

r = Cx + Du 

where x is the state vector, u is the input (or control) vector, r is the output vector, and A, 
B, C, and D are the matrix coefficients. The equations for the dynamics of the state vari- 
ables, and for the outputs of the aircraft dynamics, sensors, logic, and actuators are given 
in the following sections. The equations are then combined with equations for aerodynamics 
of the aircraft, and the total system matrix equations are formulated using the compatibility 
relationships among the dynamic systems. 

Aircraft Dynamics 

The FCAP aircraft dynamics equations for N^degrees of freedom (six rigid-body and 
(N-6) flexible-body degrees of freedom) are restricted to small perturbations which reduce 
the equations to linear form. This is a reasonable approximation for ACT studies (e.g., refs. 
4 and 5). The aircraft equation of motion expressed in state vector form is (ref. 1): 

" ^RR^R '^RD’^D ^R "R 

where 

X|^ = state vector of the displacement variables = lx,y,z,#,0,Y,qp-, . . . q|^J^ 

X|^ = state vector of the rate variables = [u, v,w,p,q,r,q^, . . . q^^ 
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Also, M|^ = generalized mass matrix, 
gravity-force matrix, = generalized aerodynamic forces in uniform flow, and u|^ = forces 
due to turbulence. 

The output of the aircraft dynamics at the sensor locations may be expressed linearly in 
terms of displacements, Xp^, rates, x^, and accelerations, x^, and, therefore, it is possible 
to write UK K 

■'D " ^SD ^ *^SR ^ *^SR ^R 

where the coefficient matrices Ucr>/ Ucn/ and UJ-d are functions of the types of sensors and 
their location. 


Arr - Coriolis force/damping matrix, Aj^p = stiffness/ 


Control System Dynamics 


The control system is defined as consisting of sensors, control logic, and actuators for 
FCAP. Classically, control system dynamics are expressed in the form of a transfer function 
which can be redefined in the state vector form of Equation (1). In the following sections, 
let x^, X|^ and x^ be the state vectors for the sensors, logic, and actuators, respectively. 

Sensors 


The state vector equations for sensor dynamics are given by 

^ss 


and 


'‘S 

where u^ is the input to the sensor system from the aircraft. 
Control Logic 


(4) 

(5) 


The state vector equations for the control logic dynamics are expressed 

x^ - A|^|^ X|^ + B|^ (U|^ + Uj^) (6) 

and 

ri = X|_ + D|_ (U|_ + u[) (7) 

where U|^ is the input to the control logic from the sensors, and ul is pilot and guidance 
system commands (see fig. 1). ^ 
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Actuators 


The equations for the dynamics of the actuators in FCAP are given by 

“ '^AA '‘A + ®A "A (®) 

and 

= <=A ''A W 

where u^ is the input to the actuators. 

Note that a term of the type u^ is absent in Equation (9). This implies that in the 
transfer function of the actuators the degree of the numerator is lower than the degree of 
the denominator. This yields considerable advantage in expressing the low-frequency- 
aerodynamics closed-loop system. 


Aerodynamics 


The potential aerodynamic method developed in references 6 to 9 provides a unified 
approach for both steady and unsteady subsonic and supersonic aerodynamics around complex, 
three-dimensional configurations. The subsonic portion of this method is incorporated into 
FCAP. The aerodynamic method of reference 6 is compatible with FCAP in that the gener- 
alized aerodynamic forces are proportional to the aircraft dynamic generalized coordinates, 
Xp, and rates, Xn, and to actuator (i.e., control surface) deflections, r^. In the time 
domain, the unsteady aerodynamic forces are expressed as 

"^RD ^ *^RA ''A^ 


where U[^[^, and U[^a are operators corresponding to frequency-dependent matrices 

Orr, andUi^^ usually known as aerodynamic-influence-coefficient matrices. 

For system stability and performance analyses, low-frequency aerodynamics is often 
adequate. Therefore, in the rest of this paper only low-frequency aerodynamics is con- 
sidered. In this case, the equations for the aerodynamics become linear with constant 
coefficients. The low-frequency aerodynamics equations are expressed in the time domain as 


j(0), 

^RR ''R 


Up = q(UX';^x„ + U)."r4 + 


j(0) 

^RD 


^RA '^A^ 


+ q(u^RRXR + *^RD^D ■'a^ 

where, for example, and 1)1^ are the first two terms of the MacLaurin-Taylor series 
of U|^R* All of these coefficienrs are frequency-independent. 
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ACT System 


The state vector equations for each of the ACT subsystems are given by Equations (2) to 
(9). Using the low-frequency aerodynamics given by Equation (11) and recognizing from 
Figure 1 that the input for each subsystem is the output of the previous subsystem, the ACT 
system can be cast in the form of Equation (1) as 


X = Ax + Bu 


where 


X = [x 


D ^A^ 


(12) 

(13) 


and 


A = F"^ A^ 
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The input u to the ACT system depends upon the gust forces, u|^, and the pilot and guidance 
system commands, u|^. 

FCAP PERFORMANCE ANALYSIS 

Performance analysis routines are available In FCAP to compute frequency response, 
transfer function poles and zeros, mean-square response to random inputs and time response. 
The theoretical basis for each of these techniques is v/ell founded in the literature; there- 
fore, the following discussions will emphasize the nature of the technique as applied to 
FCAP equations. 


Frequency Response 


Since the low-frequency-aerodynamics ACT system dynamics are expressed in state 
vector format in terms of the (constant) A, B, C, D matrices, the frequency response of the 
|<th output to the ^^'h input can be obtained (in the s-plane with zero initial conditions) from 

(s I - A) X = u^ (19) 


and 


\ = ^ < 20 ) 

Solving Equation (19) for x and substituting into Equation (20) yields 

fik/') (21) 

Classical frequency response is obtained for the special case where s = iui by computing the 
amplitude and phase from Equation (21) for a range of frequencies. 

Poles and Zeros 


Poles and zeros ate evaluated in FCAP using a different form of Hi than that given in 
Equation (21). Note that Equation (21) may be rewritten as 




Det 


D, 


k& ] '^k 


B 


I 

I i 


si - A 



( 22 ) 


The poles of H|^£ are the zeros of the denominator, i.e., the eigenvalues of the matrix A. 
The zeros of are the zeros of the numerator. The procedure to obtain them is given in 
references 1 and 10. 
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Mean-Square Response to Random Inputs 

The response of a flexible aircraft to stationary random inputs may be described in the 
frequency domain in terms of the spectral density matrix R(u3). The definition used for the 
outputs are also valid for the inputs and the state variables. 

Whej;e the inputs are represented as zero-mean white noise with constant spectral density 
matrix, U, the covariance matrix of the outputs R(0) is given (for D = 0) by 

R(0) = CX(0)C^ (23) 

where the covariance matrix of the state variables, X(0), is determined from the linear 
matrix equation (ref. 11, pp. 330-332) 

AX(0) + X(0) + BU(0) = 0 (24) 

Time Response 

The time response of a closed-loop system with initial conditions, x(0), and an arbitrary 
input function can be determined in FCAP using a fourth-order Runge-Kutta numerical inte- 
gration algorithm. For the special case where the input for t > 0 has a rational Laplace 
transform (i.e., the input can be described as a transfer function), the solution technique 
described in reference 12 is used. 

CONCLUDING REMARKS 

The equations of motion for a flexible aircraft have been presented in matrix format and 
incorporated into a digital computer program FCAP. The objective in the development of 
FCAP was to model realistically those factors that significantly affect the stability and re- 
sponse of a flexible, ACT-configured vehicle. It should be noted, however, that FCAP is 
intended primarily for use in the analysis of the performance of an ACT system, rather than 
in the synthesis of the control system. 

The small-perturbation equations of motion for the ACT system were obtained in state- 
vector format and were completed by the addition of aerodynamics of arbitrary, complex 
aircraft configurations. Both fully unsteady and low-frequency aerodynamic equations were 
presen^ted; however, for performance analyses, low-frequency aerodynamics is usually 
adequate and, therefore, the ACT system equations were presented for low-frequency aero- 
dynamics only. Progrom FCAP, however, allows the analysis of complete aircraft configu- 
rations, including control system, with either low-frequency or unsteady aerodynamics. 

The analysis of ACT system performance in FCAP has also been presented. In particular, 
the application of classical frequency response, poles and zeros, mean-square response, and 
time response in FCAP in state-vector format has been discussed. 

Program FCAP provides a computerized method of integrating multiple systems into a 
matrix format, and then provides the means for obtaining desired solutions through classical 
analysis techniques. The program is currently in the final stages of checkout and has been 
used to solve textbook examples of control system problems. The program thus far has proven 
to be simple to use and requires a minimum of input. 
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Figure 1. ACT System. 










